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Abstract. Phase transitions of atmospheric water play a 
ubiquitous role in the Earth's climate system, but their direct 
impact on atmospheric dynamics has escaped wide attention. 
Here we examine and advance a theory as to how conden- 
sation influences atmospheric pressure through the mass re- 
moval of water from the gas phase with a simultaneous ac- 
count of the latent heat release. Building from the fundamen- 
tal physical principles we show that condensation is associ- 
ated with a decline in air pressure in the lower atmosphere. 
This decline occurs up to a certain height, which ranges from 
3 to 4 km for surface temperatures from 10 to 30 °C. We 
then estimate the horizontal pressure differences associated 
with water vapor condensation and find that these are com- 
parable in magnitude with the pressure differences driving 
observed circulation patterns. The water vapor delivered to 
the atmosphere via evaporation represents a store of poten- 
tial energy available to accelerate air and thus drive winds. 
Our estimates suggest that the global mean power at which 
this potential energy is released by condensation is around 
one per cent of the global solar power - this is similar to the 
known stationary dissipative power of general atmospheric 
circulation. We conclude that condensation and evaporation 
merit attention as major, if previously overlooked, factors in 
driving atmospheric dynamics. 



1 Introduction 

Phase transitions of water are among the major physical pro- 
cesses that shape the Earth's climate. But such processes 
have not been well characterized. This shortfall is recognized 
both as a challenge and a prospect for advancing our un- 
derstanding of atmospheric circulation (e.g. 



Lorenz 1983 



| Schneider] |2006[ ). In A History of Prevailing Ideas about the 
General Circulation of the Afmo^p/zere |Lorenz| ([T983 ) wrote: 
"We may therefore pause and ask ourselves whether this step will 
be completed in the manner of the last three. Will the next decade 
see new observational data that will disprove our present ideas? It 
would be difficult to show that this cannot happen. 

Our current knowledge of the role of the various phases of wa- 
ter in the atmosphere is somewhat incomplete: eventually it must 
encompass both thermodynamic and radiational effects. We do not 
fully understand the interconnections between the tropics, which 
contain the bulk of water, and the remaining latitudes. . . . Perhaps 
near the end of the 20th century we shall suddenly discover that we 
are beginning the fifth step." 
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that local pressure is reduced by precipitation and increased 
by evaporation. |Qiu et aL] ([1993 ) noted that "the mass de- 
pletion due to precipitation tends to reduce surface pressure, 
which may in turn enhance the low-level moisture conver- 
gence and give a positive feedback to precipitation". [Van 
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den Dool and Saha ( 1993) labeled the effect as a physically 
distinct "water vapor forcing". |LackmantT and Yab lonsky] 
( |2004[ ) investigated the precipitation mass sink for the case 
of Hurricane Lili (2002) and made an important observation 
that "the amount of atmospheric mass removed via precip- 
itation exceeded that needed to explain the model sea level 
pressure decrease". 

Although the pressure changes associated with evapora- 
tion and condensation have received some attention, the in- 
vestigations have been limited: the effects remain poorly 
characterized in both theory and observations. Previous in- 
vestigations focused on temporal pressure changes not spa- 
tial gradients. Even some very basic relationships remain 
subject to confusion. For example, there is doubt as to 
whether condensation leads to reduced or to increased atmo- 
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spheric pressure (Poschl 2009 p. S12436). Opining that the 
status of the issue in the meteorological literature is unclear, 



Haynes (2009) suggested that to justify the claim of pres- 



sure reduction one would need to show that "the standard 
approaches (e.g., set out in textbooks such as "Thermody- 
namics of Atmospheres and Oceans" by ICurry and Webster 



1999)) imply a drop in pressure associated with condensa- 



tion". 

Here we aim to clarify and describe, building from basic 
and established physical principles, the pressure changes as- 
sociated with condensation. We will argue that atmospheric 
water vapor represents a store of potential energy that be- 
comes available to accelerate air as the vapor condenses. 
Evaporation, driven by the sun, continuously replenishes the 
store of this energy in the atmosphere. 

The paper is structured as follows. In Section 2 we an- 
alyze the process of adiabatic condensation to show that it 
is always accompanied by a local decrease of air pressure. 
In Section 3 we evaluate the effects of water mass removal 
and lapse rate change upon condensation in a vertical air col- 
umn in approximate hydrostatic equilibrium. In Section 4 we 
estimate the horizontal pressure gradients induced by water 
vapor condensation to show that these are sufficient enough 
to drive the major circulation patterns on Earth (Section 4.1). 
We examine why the key relationships have remained un- 
known until recently (Section 4.2). We evaluate the mean 
global power available from condensation to drive the gen- 
eral atmospheric circulation (Secton4.3). Finally, we discuss 
the interplay between evaporation and condensation and the 
essentially different implications of their physics for atmo- 
spheric dynamics (Section 4.4). In the concluding section 
we discuss the importance of condensation as compared to 
differential heating as the major driver of atmospheric cir- 
culation. Our theoretical investigations strongly suggest that 
the phase transitions of water vapor play a far greater role in 
driving atmospheric dynamics than is currently recognized. 

2 Condensation in a local air volume 

2.1 Adiabatic condensation 

We will first show that adiabatic condensation is always ac- 
companied by a decrease of air pressure in the local vol- 
ume where it occurs. The first law of thermodynamics 
for moist air saturated with water vapor reads ( |Gill| |1982| 
|Curry and Webster] [1999] ) 

dQ = c v dT + pdV + Ld'y, (1) 

Pv ^ 1 dj dp v dp 

7 = — < 1, — = . (2) 

P 7 Pv P 

Here p v is partial pressure of saturated water vapor, p is air 
pressure, T is absolute temperature, Q (J mol -1 ) is molar 
heat, V (m 3 mol -1 ) is molar volume, L « 45 kJ mol -1 
is the molar heat of vaporization, cy — |i? is molar heat 



capacity of air at constant volume (J mol -1 K _1 ), R = 
8.3 J mol -1 K _1 is the universal gas constant. In processes 
not involving phase transitions the third term in ([TJ is zero. 
In such processes partial pressure p v changes proportionally 
to air pressure p, so that function 7 |2} does not change. The 
small value of 7 < 0.1 under terrestrial conditions allows us 
to neglect the influence made by the heat capacity of liquid 
water in Eq. ([TJ. 

The partial pressure of saturated water vapor obeys the 
Clausius-Clapeyron equation: 



dp v _ dT 
Pv T ' 



L 

RY 



(3) 



Pv (T) = p v0 exp(£ - , ( 4 ) 

where p„o an d £0 correspond to some reference temperature 
T Q . Below we use T = 303 K and p vQ = 42 hPa ( |Bolton| 
1980) and neglect the dependence of L on temperature. 

We will also use the ideal gas law as the equation of state 
for atmospheric air: 

pV = RT, (5) 
dp dV dT 

— + ^ = ^r- (6) 
p V 1 

Using Eq. (|6]l the first two terms in Eq. ([TJl can be written 
in the following form 



cydT + pdV 



RT 
0.29 



dT 

Y 



dp 

- V — 
P 

c v + R. 



(7) 



Writing dj in ([TJ with use of Q and ((3J as 
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dp 
P 



(8) 



and using the definition of £ <|3j we arrive at the following 
form for the first law of thermodynamics ([TJ: 



dQ = 



RT 



T p 



(9) 



In adiabatic processes dQ = 0, and the expression in 
braces in Q turns to zero, which implies: 



T p 



1+7C 

1 + vie 



(10) 



Note that p, 7 and £ are all dimensionless; 7 and £ are vari- 
ables and /j, is a constant, tp(0, 0) = fi. This is a general 
dependence of temperature on pressure in an adiabatic atmo- 
spheric process that involves phase transitions of water vapor 
(evaporation or condensation), i.e. change of 7. At the same 
time 7 itself is a function of temperature as determined by 
Eq. (§: 



e?7 

7 



.dT 



dp 
P 



dT(<p-l 



T 



dp 



(11) 
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One can see from Eqs. ( [T0| and (Hi that the adiabatic phase 
transitions of water vapor are fully described by the relative 
change of either pressure dp/p or temperature dT/T. For the 
temperature range relevant for Earth we have £ = Lj RT w 
18 so that 

1-4.3. (12) 



dp < 0). We write the first law of thermodynamics (|9|l and 
Eq. (fTT 



dS 



R 1 + /n£ 2 



dT dp dT _ 1 /d7 
p ' T £ V 7 



dp 
P 



Excluding dT/T from Eqs. ( 17 1 we obtain 



Noting that p, 7, £ are all positive, from ( 10 1, ( 1 1 1 and ( 12 1 
we obtain 

1 + 7f f p - 1 

&>-l = frt , ...^ -1= , >0- (13) 



dp 



- 1) 



7 



(17) 



(18) 



1 + Pit 2 



1 + ^ 2 



p 71+ ^ i? 

The term in round brackets in Eq. (jT8j is positive, see ( pj ) 
the multiplier at dS is also positive. Therefore, when con- 
densation occurs, i.e., when g?7/7 < 0, and dS > 0, the left 



Condensation of water vapor corresponds to a decrease of hand side G f Eq. tpl is negative. This means that dp/p < 0, 



7, dj < 0. It follows unambiguously from Eqs. (Hi and ( 13 1 
that if dj is negative, then dp is negative too. This proves that 
water vapor condensation in any adiabatic process is neces- 
sarily accompanied by reduced air pressure. 

2.2 Adiabatic condensation cannot occur at constant vol- 
ume 

Our previous result refutes the proposition that adiabatic con- 
densation can lead to a pressure rise due to the release of 



latent heat ( |Poschl| |2009] p. S 12436). Next, we show that 
while such a pressure rise was implied by calculations as- 
suming adiabatic condensation at constant volume, in fact 
such a process is prohibited by the laws of thermodynamics 
and thus cannot occur. 

Using |6} and ( fTO] ) we can express the relative change of 
molar volume dV/V in terms of dy/y: 



dV 



1 — ip dy 

<p£ - 1 7 ' 



Putting dV = in ( p| we obtain 

(1 - tp)dy 



(14) 



(15) 



The denominator in (|T5]i is greater than zero, see Eq. ( 13 I. 
In the numerator we note from the definition of ip ( fT0] i that 



l-<p = 



- 27 

7+2 7 f 2 



_5_ 

27 



The expression in square 



brackets lacks real roots: 



7 < 1. (16) 



In consequence, Eq. ( 15 1 has a single solution dj = 0. This 
proves that condensation cannot occur adiabatically at con- 
stant volume. 

2.3 Non-adiabatic condensation 

To conclude this section, we show that for any process where 
entropy increases, dS — dQ /T > 0, water vapor condensa- 
tion (g?7 < 0) is accompanied by drop of air pressure (i.e., 



i.e., air pressure decreases. 

Condensation can be accompanied by a pressure increase 
only if dS < 0. This requires that work is performed 
on the gas such as occurs if it is isothermally compressed. 
(We note too, that if pure saturated water vapor is isother- 
mally compressed condensation occurs, but the Clausius- 
Clapeyron equation ([3]l shows that the vapor pressure re- 
mains unchanged being purely a function of temperature.) 



3 Adiabatic condensation in the gravitational field 

3.1 Difference in the effects of mass removal and tempera- 
ture change on gas pressure in hydrostatic equilibrium 

We have shown that adiabatic condensation in any local vol- 
ume is always accompanied by a drop of air pressure. We 
will now explore the consequences of condensation for the 
vertical air column. 

Most circulation patterns on Earth are much wider than 
high, with the ratio height/length being in the order of 10~ 2 
for hurricanes and down to 10 -3 and below in larger regional 
circulations. As a consequence of mass balance, vertical ve- 
locity is smaller than horizontal velocities by a similar ra- 
tio. Accordingly, the local pressure imbalances and result- 
ing atmospheric accelerations are much smaller in the verti- 
cal orientation than in the horizontal plane, the result being 
an atmosphere in approximate hydrostatic equilibrium ( |Gill| 
|1982| l. Air pressure then conforms to the equation 

dp 



dz 



= P9, P(0) = Ps 



g j p(z)dz. 
'0 



(19) 



Applying the ideal gas equation of state <|3j we have from 

%. = p h =xr m (20) 

dz h Mg 



This solves as 



p(z) = p s exp <^ - 



dz 1 

W) 



(21) 



Here M is air molar mass (kg mol x ), which, as well as 
temperature T{z), in the general case also depends on z. 
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The value of p s (p~9]>, air pressure at the surface, appears 



as the constant of integration after Eq. ( 19 1 is integrated over 
z. It is equal to the weight of air molecules in the atmo- 
spheric column. It is important to bear in mind that p s does 
not depend on temperature, but only on the amount of gas 
molecules in the column. It follows from this observation 
that any reduction of gas content in the column reduces sur- 
face pressure. 

Latent heat released when water condenses means that 
more energy has to be removed from a given volume of satu- 
rated air for a similar decline in temperature when compared 
to dry air. This is why the moist adiabatic lapse rate is smaller 
than the dry adiabatic lapse rate. Accordingly, given one and 
the same surface temperature T s in a column with rising air, 
the temperature at some distance above the surface will be 
on average higher in a column of moist saturated air than in 
a dry one. 

However, this does not mean that at a given height air pres- 
sure in the warmer column is greater than air pressure in the 
colder column (cf. Meesters et al. 2009; Makarieva and Gor- 
shkov 2009c|, because air pressure p(z) ( 21 1 depends on two 



parameters, temperature T(z) and surface air pressure (i.e., 
the total amount of air in the column). If the total amount of 
air in the warmer column is smaller than in the colder col- 
umn, air pressure in the surface layer will be lower in the 
warmer column despite its higher temperature. 

In the following we estimate the cumulative effect of gas 
content and lapse rate changes upon condensation. 

3.2 Moist adiabatic temperature profile 

Relative water vapor content d2|) and temperature T depend 



on height z. From Eqs. ( [TO i 

dT 



T 



and ( 20 1 we have 



1 dry £<y? — 1 
7 dz h 



Cm- i i 

1 + my£ 2 h' 



(22) 



(23) 



Eq. ( 22 1 represents the well-known formula for moist adia- 
batic gradient as given in Glickman (2000) for small 7 < 0.1. 
At 7 = we have <p(7, £) = \i and Td = M^g/cp = 
9.8 K km -1 , which is the dry adiabatic lapse rate that is inde- 
pendent of height z, = 29 g mol -1 . For moist saturated 
air the change of temperature T and relative partial pressure 
7 of water vapor with height is determined by the system of 
differential equations <\22\ , ( |23"j ). 

Differentiating both parts of Clapeyron-Clausius equation 
([3]l over z we have, see ( |22] ): 



dp v 
dz 



Pv_ 



h, 



RT 2 



LT 



T 



Pv(z) = p vs exp 



dz' 



Pvs 



h 

Pv(0). 



(24) 



The value of h v represents a fundamental scale height for the 
vertical distribution of saturated water vapor. At T s = 300 K 
this height h v is approximately 4.5 km. 

Differentiating both parts of Eq. Q over z with use of ( f20| > 
and (24 1 and noticing that h v = h/(£<p) we have 



1 d-y 

7 dz 



1 dp v 

p v dz 



I dp 
p p 



1 

K 



1 

h 



h v h 

h — h v 
(25) 



This equation is equivalent to Eq. (23 1 when Eqs. (22 1 and 
(24 1 are taken into account. Height h 1 represents the vertical 



scale of the condensation process. Height scales h v (24i and 



h~ (25 1 depend on 95(7, £) (22 1 and, consequently, on 7. At 



T s = 300 K height /i 7 w 9 km, in close proximity to the 



water vapor scale height described by Mapes (2001 ). 



3.3 Pressure profiles in moist versus dry air columns 

We start by considering two static vertically isothermal at- 
mospheric columns of unit area, A and B, with temperature 
T(z) — T s independent of height. Column A contains moist 
air with water vapor saturated at the surface, column B con- 
tains dry air only. Surface temperatures and surface pressures 
in the two columns are equal. In static air Eq. ( p"9| ) is exact 
and applies to each component of the gas mixture as well as 
to the mixture as a whole. At equal surface pressures, the 
total air mass and air weight are therefore the same in both 
columns. Water vapor in column A is saturated at the sur- 
face (i.e., at z = 0) but non-saturated above it (at z > 0). 
The saturated partial pressure of water vapor at the surface 
Pv(T s ) Q is determined by surface temperature and, as it is 
in hydrostatic equilibrium, equals the weight of water vapor 
in the static column. 

We now introduce a non-zero lapse rate to both columns: 



the moist adiabatic V (22i to column A and the dry adia- 
batic in column B. (Now the columns cannot be static: the 
adiabatic lapse rates are maintained by the adiabatically as- 
cending air.) Due to the decrease of temperature with height, 
some water vapor in column A undergoes condensation. Wa- 
ter vapor becomes saturated everywhere in the column (i.e., 
at z > 0), with pressure p v (z) following Eq. (24 1 and density 
p v = p v M v /(RT) =p v / (gh n ) following 



p v (z) = p v (T s ) 

Pv(T s ) = , K 

gh n (T s ) 



exp 



h n (z) 
_ RT(z) 
= M v g 



dz' 

h v (z') j ' 
T(z) - T s - Tz. 

(26) 



Here h n (z) is the scale height of the hydrostatic distribution 
of water vapor in the isothermal atmosphere with T s = T(z). 

The change in pressure 5p s in column A due to water vapor 
condensation is equal to the difference between the initial 
weight of water vapor p v (T s ) and the weight of saturated 
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Fig. 1. (a): scale height of saturated water vapor h v (z) {24}, hydrostatic scale height of water vapor h n (z) 1 26 i, and scale height of moist 



air h(z) q20p in the column with moist adiabatic lapse rate {22} for three values of surface temperature T 3 ; (b): condensation-induced drop 
of air pressure at the surface {27} as dependent on surface temperature T s ; (c): pressure difference versus altitude z between atmospheric 
columns A and B with moist and dry adiabatic lapse rates, Eqs. {30}, l |31[ >, respectively, for three values of surface temperature T s . Height z c 
at which pa(z c ) — Pb(z c ) = is 2.9, 3.4 and 4.1 km for 283, 293 and 303 K, respectively. Due to condensation, at altitudes below z c the 
air pressure is lower in column A despite it being warmer than column B. 



water vapor: 



gravitational acceleration (giving weight per unit area): 



Sp s = Pv(T s ) -g J p v (z)dz < p v (T s ) - p v {T s )gh v (T s ) 

(27) 



The inequality in Eq. ( 27 1 represents a conservative estimate 
of Sp s due to the approximation h v (z) — h v (T s ) made while 



integrating p v (z) (26i. As far as h v (z) declines with height 



more rapidly than h n {z), Fig. la, the exact magnitude of this 
integral is smaller, while the value of Sp s is larger. The phys- 



ical meaning of estimate ( 27 1 consists in the fact that the drop 



of temperature with height compresses the water vapor distri- 
bution h ns /h vs -fold compared to the hydrostatic distribution 
( [Makarieva and Gorshkov| [20071 |2009a) . 

The value of Sp s (27 1 was calculated as the difference 
between the weight per unit surface area of vapor in the 
isothermal hydrostatic column and the weight of water va- 
por that condensed when a moist adiabatic lapse rate was 
applied. This derivation can also be understood in terms of 
the variable conventionally called the adiabatic liquid wa- 
ter content (e.g., Curry and Webster 1999 Eq. 6.41). We 



can represent the total mixing ratio of moisture (by mass) as 
q t = q v + qi = p v /p + pi/p, where p v is the mass of va- 
por gj c 1 and pi is the mass of liquid water per unit air 
volume; qt <C 1. The total adiabatic liquid water content 
in the column equals the integral of qip over z at constant 
It, 1W = QtP - q v R = qtP - Pv The value of 6p s {27} is 
equal to this integral (mass per unit area) multiplied by the 



5p s 



g I qipdz 
'o 



q t pdz 



Pvdz . (28) 



The first integral in the right-hand part of this equation gives 
the mass of vapor in the considered atmospheric column if 
water vapor were a non-condensable gas, q v = q t = const. 
This term is analagous to the first term, p v (T s ), in the right- 



hand side of Eq. ( 27 1, where a static isothermal column was 
considered. The second term is identical to the second term, 



g J Q p v dz, in Eq. (27 1. 

Using the definition of h v (T s ) {24}, h n (T s ) {26} and re- 
calling that M v /M d = 0.62 and p v (T s ) = j s p s , see {4i, 
we obtain the following expression for the Sp s estimate ( 27 1, 
Fig. lb: 



Sp s 

Ps 



Is 1 - 0.62 



1 + IsPg 



(29) 



" fj£s + isp£, 2 s , 

Note that 5p s /p s is proportional to ^ s and increases expo- 
nentially with the rise of temperature. 

After an approximate hydrostatic equilibrium is estab- 
lished, the vertical pressure profiles for columns A and B 
become, cf. {2T}: 

Sp / [* dz' \ RT 

Ps I Jo h A (z')) Mg 

(30) 



Pb{z) = p s exp 
Here M(z) = M d {\ 



dz 1 



h B {z') 



h B = 



RT d 



(31) 



7) 



M vT , l = 
T{z) obey Eqs. {22} and {23}, T d (z) = T s 



M d g 
Pv{z)/pa{z) and 
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In Fig. lc the difference pa (z) — pb (z) is plotted for three 
surface temperatures, T s = 10°, 20° and 30° C. In all three 
cases condensation has resulted in a lower air pressure in col- 
umn A compared to column B everywhere below z c 2.9, 
3.4 and 4.1 km, respectively. It is only above that height that 
the difference in lapse rates makes pressure in the moist col- 
umn higher than in the dry column. 



4 Relevance of the condensation-induced pressure 
changes for atmospheric processes 

4.1 Horizontal pressure gradients associated with vapor 
condensation 

We have shown that condensation of water vapor produces a 
drop of air pressure in the lower atmosphere up to an altitude 
of a few kilometers, Fig. lc, in a moist saturated hydrostat- 
ically adjusted column. In the dynamic atmospheric context 
the vapor condenses and latent heat is released during the 
ascent of moist air. The vertical displacement of air is in- 
evitably accompanied by its horizontal displacement. This 
translates much of the condensation-induced pressure differ- 
ence to a horizontal pressure gradient. Indeed, as the up- 
welling air loses its water vapor, the surface pressure dimin- 
ishes via hydrostatic adjustment producing a surface gradi- 
ent of total air pressure between the areas of ascent and de- 
scent. The resulting horizontal pressure gradient is propor- 
tional to the the ratio of vertical to horizontal velocity w/u 
( [Makarieva and Gorshkov||2009b"l l. 

We will illustrate this point regarding the magnitude of the 
resulting atmospheric pressure gradient for the case of a sta- 
tionary axis-symmetric circulation developing above a hor- 
izontally isothermal oceanic surface. In cylindrical coordi- 
nates the continuity equation for the mixture of condensable 
(vapor) and non-condensable (dry air) gases can be written 
as 

ld(N d ur) d(N d w) 



S{r, z) 



w 



r dr 
1 d(N v ur) 
r dr 
^>N V 



dz 

d{N v w) 



dz 



N 



ON 

dz 



0; 

S(r, z) 
$7 



(32) 



(33) 



u>7ir> N = N V +N d . 
dz 

(34) 

Here N d and N v are molar densities of dry air and water 
vapor, respectively; 7 = N v /N, see pi, r is the distance 
from the center of the area where condensation takes place, 
S(r, z) is the sink term describing the non-conservation of 
the condensable component (water vapor). Saturated pres- 
sure of water vapor depends on temperature alone. Assum- 
ing that vapor is saturated at the isothermal surface we have 
8N v /dr = 0, so N v only depends on z. (Note that this 
condition necessitates either that there is an influx of water 
vapor via evaporation from the surface (if the circulation pat- 
tern is immobile), or that the pressure field moves as vapor is 



locally depleted. The second case occurs in compact circu- 
lation patterns like hurricanes and tornadoes^]) As the air as- 
cends with vertical velocity w, vapor molar density decreases 
due to condensation and due to the expansion of the gas along 
the vertical gradient of decreasing pressure. The latter ef- 
fect equally influences all gases, both condensable and non- 
condensable. Therefore, the volume-specific rate S(r, z) at 
which vapor molecules are locally removed from the gaseous 
phase is equal to w(dN v /dz - (N v /N)dN/dz), see (jT), ®. 
The second term describes the expansion of vapor at a con- 
stant mixing ratio which would have occurred if vapor were 
non-condensable as the other gases. (If vapor did not con- 
dense, its density would decrease with height as a constant 
proportion of the total molar density of moist air as with any 
other atmospheric gas.) 



The mass of dry air is conserved, Eq. ( 32 1. Using this fact. 



Eq. ( 34 1 and dN v /dr = one can see that 



N 



1 d(ur) 
r dr 



dw 
dz 



dN 
dz 



0. 



(35) 



Now expressing dN/dr — dN d /dr + dN v /dr from 



Eqs. ( 32 1 and (B3JI with use of Eq. p5f we obtain 



dN _ w fdN v N v dN 
dr u \ dz N dz 



(36) 



Using the equation of state for moist air p = NRT and water 



vapor p v — N V RT we obtain from Eqs. ( 36 1 and ( 25 1: 



dp 
dr 



dp v 
dz 



Pv_dp 
p dz 



w 
u 



7P w 
hj u 



(37) 



Here velocities w and u represent vertical and radial veloci- 
ties of the ascending air flow, respectively. The ascending air 
converges towards the center of the area where condensation 
occurs. Scale height h 7 is defined in Eq. |25] l. A closely re- 
lated formula for horizontal pressure gradient can be applied 
to a linear two-dimensional air flow, with dp/dr replaced by 
dp/dx. 



Equation ( 37 1 shows that the difference between the scale 



heights h v and h ( 25 1 of the vertical pressure distributions 



for water vapor and moist air leads to the appearance of a 
horizontal pressure gradient of moist air as a whole. This 
equation contains the ratio of vertical to horizontal velocity. 
Estimating this ratio it is possible to evaluate, for a given cir- 
culation, what sorts of horizontal pressure gradients are pro- 
duced by condensation and whether these gradients are large 
enough to maintain the observed velocities via the positive 
physical feedback described by Eq. ( [37) . 

For example, for Hadley cells at T = 300 K, h~ = 9 km, 
7 = 0.04 and a typical ratio of w/u ~ 10 -3 we obtain from 



'Makarieva, A. M. and Gorshkov, V. G.: Potential en- 
ergy of atmospheric water vapor and the air motions in- 
duced by water vapor condensation on different spatial scales, 
http://arxiv.org/abs/1003.5466 [physics.gen-ph], 2010. 
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Eq. ( 37 i a pressure gradient of about 0.4 Pa km -1 . On a dis- 
tance of 3000 km such a gradient would correspond to a pres- 
sure difference of 12 hPa, which is close to the upper range of 
the actually observed pressure differences in the region (e.g., 
Murp hree and Van den Dool"}|1988| Fig. 1). This estimate 



illustrates our proposal that condensation should be consid- 
ered one of the main determinants of atmospheric pressure 
gradients and, hence, air circulation. 

Similar pressure differences and gradients, also compara- 



ble in magnitude to Sp s (27 1 and dp/dr (37i are observed 
within cyclones, both tropical and extratropical, and persis- 
tent atmospheric patterns in the low latitudes (Holland 1980; 
|Zhou and Laul[T9l^|B7u"mmer et al.|[2000l|NicholsonH200o[ 



Simmonds et al. 20081. For example, the mean depth of 



Arctic cyclones, 5 hPa (Simmonds et al. 20081, is about 
ten times smaller than the mean depth of a typical tropi- 
cal cyclone (Holland"! |1980| l. This pattern agrees well with 
the Clausius-Clapeyron dependence of Sp s , Fig. lb, which 
would predict an 8 to 16-fold decrease with mean oceanic 
temperature dropping by 30-40 degrees Celsius. The exact 
magnitude of pressure gradient and the resulting velocities 
will depend on the horizontal size of the circulation pattern, 
the magnitude of friction and degre e of the radial symmetry 
( |Makarieva and Gorshkov 2009a bp 1 . 



4.2 Regarding previous oversight of the effect 

For many readers a major barrier to acceptance of our propo- 
sitions may be to understand how such a fundamental physi- 
cal mechanism has been overlooked until now. Why has this 
theory come to light only now in what is widely regarded as 
a mature field? We can offer a few thoughts based on our 
readings and discussions with colleagues. 

The condensation-induced pressure gradients that we have 
been examining are associated with density gradients that 
have been conventionally considered as minor and thus ig- 
nored in the continuity equation (e.g., Sabato, 2008). For 
example, a typical Ap = 50 hPa pressure difference ob- 
served along the horizontally isothermal surface between the 
outer environment and the hurricane center (e.g., |Holland| 
1980 ) is associated with a density difference of only around 



5%. This density difference can be safely neglected when 
estimating the resulting air velocity u from the known pres- 
sure differences Ap. Here the basic scale relation is given 
by Bernoulli's equation, pu 2 /2 = Ap. The point is that 
a 5% change in p does not significantly impact the magni- 
tude of the estimated air velocity at a given Ap. But, as we 
have shown in the previous section, for the determination of 



the pressure gradient ( 37 1 the density difference and gradient 
< [36] > are key. 

Considering the equation of state |5]l for the horizontally 
isothermal surface we have p = Cp, where C = RT/M = 
const. Irrespective of why the considered pressure differ- 
ence arises, from Bernoulli's equation we know that u 2 — 
2 Ap/p = 2CAp/p, Ap = p — p. Thus, if one puts 



Ap/ p = Ap/p equal to zero, no velocity forms and there is 
no circulation. Indeed, we have u 2 — 2 Ap/p — 2CAp/p = 
2C(Ap/ po)(l + Ap/ po + ...). As one can see, discarding 
Ap compared to p does indeed correspond to discarding the 
higher order term of the smallness parameter Ap/ p. But with 
respect to the pressure gradient, the main effect is propor- 
tional to the smallness parameter Ap/ po itself. If the latter 
is assumed to be zero, the effect is overlooked. We suggest 
that this dual aspect of the magnitude of condensation-related 
density changes has not been recognized and this has con- 
tributed to the neglect of condensation-associated pressure 
gradients in the Earth's atmosphere. 

Furthermore, the consideration of air flow associated with 
phase transitions of water vapor has been conventionally re- 
duced to the consideration of the net fluxes of matter. Sup- 
pose we have a linear circulation pattern divided into the as- 
cending and descending parts, with similar evaporation rates 
E (kg H2O m~ 2 s _1 ) in both regions. In the region of ascent 
the water vapor precipitates at a rate P. This creates a mass 
sink E — P, which has to be balanced by water vapor import 
from the region of descent. Approximating the two regions 
as boxes of height h, length I and width d, the horizontal ve- 
locity u t associated with this mass transport can be estimated 
from the mass balance equation 



(P — E)ld = u t phd, u t 



(P-E) l 
p h 



(38) 



Equation ( |38| l says that the depletion of air mass in the region 
of ascent at a total rate of (P~E)ld is compensated for by the 
horizontal air influx from the region of descent that goes with 
velocity u t via vertical cross-section of area hd. For typical 
values in the tropics with P — E ~ 5 mm day -1 = 5.8 x 
10" 5 kg H 2 m~ 2 s" 1 and l/h ~ 2 X 10 3 we obtain u t ~ 
1 cm s _1 . For regions where precipitation and evaporation 
are smaller, th e value of u t will be smaller too. For example, 
Lorenz ( 1967 p. 51) estimated u t to be ~ 0.3 cm s _1 for the 
air flow across latitude 40°S. 

With p ss pd the value of u t can be understood as the 
mass-weighted horizontal velocity of the dry air + water va- 
por mixture, which is the so-called barycentric velocity, see, 



e.g., \ Wacker and Herbertj |2003) | Wacker etal) |2006) . There 
is no net flux of dry air between the regions of ascent and de- 
scent, but there is a net flux of water vapor from the region of 
descent to the region of ascent. This leads to the appearance 
of a non-zero horizontal velocity u t directed towards the re- 
gion of ascent. Similarly, vertical barycentric velocity at the 
surface is Wt ~ (E — P)/p (Wacker and Herbert 2003 1, 



which reflects the fact that there is no net flux of dry air via 
the Earth's surface, while water vapor is added via evapo- 
ration or removed through precipitation. The absolute mag- 
nitude of vertical barycentric velocity wt for the calculated 
tropical means is vanishingly small, wt ~ 0.05 mm s _1 . 

We speculate that the low magnitude of barycentric ve- 
locities has contributed to the judgement that water's phase 
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transitions cannot be a major driver of atmospheric dynam- 
ics. However, barycentric velocities should not be confused 



(e.g., Meesters et al. 2009 1 with the actual air velocities. Un- 
like the former, the latter cannot be estimated without consid- 



ering atmospheric pressure gradients (Makarieva and Gor- 
2009c[ ). For example, in the absence of friction, the 



shkov 



maximum linear velocity u c that could be produced by con- 
densation in a linear circulation pattern in the tropics consti- 
tutes 

u c = y / 2Ap//>~ 45ms~ 1 > u t . (39) 

Here Ap was taken equal to 12 hPa as estimated from 
Eq. ( (37] > for Hadley cell in Section 4.1. As one can see, u c 
(39 1 is much greater than u t (38 I. As some part of poten- 
tial energy associated with the condensation-induced pres- 



sure gradient is lost to friction (Makarieva and Gorshkov 



2009a I, real air velocities observed in large-scale circulation 
are an order of magnitude smaller than u c , but still nearly 
three orders of magnitude greater than w t . 

4.3 The dynamic efficiency of the atmosphere 

We will now present another line of evidence for the impor- 
tance of condensation-induced dynamics: we shall show that 
it offers an improved understanding of the efficiency with 
which the Earth's atmosphere can convert solar energy into 
kinetic energy of air circulation. While the Earth on aver- 
age absorbs ab out I 2.4 x 1 2 W m~ 2 of solar radiation 
(Raval and Ramanathan 1989 1, only a minor part 77 ~ 10~ 2 



of this energy is converted to the kinetic power of atmo- 
spheric and oceanic movement. |Lorenz| ( |1967 p. 97) notes, 
"the determination and explanation of efficiency r\ constitute 
the fundamental observational and theoretical problems of 
atmospheric energetics". Here the condensation-induced dy- 
namics yields a relationship that is quantitative in nature and 
can be estimated directly from fundamental atmospheric pa- 
rameters. 

A pressure gradient is associated with a store of potential 
energy. The physical dimension of pressure gradient coin- 
cides with the dimension of force per unit air volume, i.e. 
1 Pa m _1 = 1 N m -3 . When an air parcel moves along the 
pressure gradient, the potential energy of the pressure field 
is converted to the kinetic energy. The dimension of pres- 
sure is identical to the dimension of energy density: 1 Pa = 



1 N mr z = 1 J irT 



As the moist air in the lower part of 



the atmospheric column rises to height /i 7 where most part 
of its water vapor condenses, the potential energy released 
amounts to approximately 5p s ( [27) 1. The potential energy re- 
leased ir v per unit mass of water vapor condensed, dimension 
J (kg H2O ) , thus becomes 



TT V (T S ) = 



Pv 



RT S 



1 - 



M v gT s 



(40) 



The global mean precipitation rate is P ~ 
10 3 kg H 2 Q m~ 2 year -1 dL'vovitchl Il979|l, global mean 



288 K and the observed mean 
6.5 Kkm" 1 dGlickmanl 12000b. 



surface temperature is T s 
tropospheric lapse rate r o 
Using these values and putting T instead of the moist adia- 
batic lapse rate T s in (|40j, we can estimate the global mean 
rate n„ = Ptt v at which the condensation-related potential 
energy is available for conversion into kinetic energy. At 
the same time we also estimate the efficiency 77 = II„/7 
of atmospheric circulation that can be generated by solar 
energy via the condensation-induced pressure gradients: 



n„ 



P7r„~3.5Wnr 2 , 77 ~ 0.01E 



(41) 



Thus, the proposed approach not only clarifies the dynamics 
of solar energy conversion to the kinetic power of air move- 
ment (solar power spent on evaporation — > condensation- 
related release of potential power — > kinetic power genera- 
tion). It does so in a quantiatively tractable manner explain- 
ing the magnitude of the dissipative power associated with 
maintaining the kinetic energy of the Earth's atmosphere. 

Our estimate of atmospheric efficiency differs fundamen- 
tally from a thermodynamic approach based on calculat- 
ing the entropy budgets under the assumption that the at- 



mosphere works as a heat engine, e. g., (|Pauluis and Held 



2002a bj, see also ( |Makarieva et al. 2010 1. The principal 
limitation of the entropy-budget approach is that while the 
upper bounds on the amount of work that could be produced 
are clarified, there is no indication regarding the degree to 
which such work is actually performed. In other words, the 
presence of an atmospheric temperature gradient is insuffi- 
cient to guarantee that mechanical work is produced. In con- 



trast, our estimate (41 1 is based on an explicit calculation of 
mechanical work derived from a defined atmospheric pres- 
sure gradient. It is, to our knowledge, the only available esti- 
mate of efficiency 77 made from the basic physical parameters 
that characterize the atmosphere. 

4.4 Evaporation and condensation 

While condensation releases the potential energy of atmo- 
spheric water vapor, evaporation, conversely, replenishes it. 
Here we briefly dwell on some salient differences between 
evaporation and condensation to complete our picture regard- 
ing how the phase transitions of water vapor generate pres- 
sure gradients. 

Evaporation requires an input of energy to overcome the 
intermolecular forces of attraction in the liquid water to free 
the water molecule to the gaseous phase, as well as to com- 
press the air. That is, work is performed against local atmo- 
spheric pressure to make space for vapor molecules that are 
being added to the atmosphere via evaporation. This work, 
associated with evaporation, is the source of potential energy 
for the condensation-induced air circulation. Upon conden- 
sation, two distinct forms of potential energy arise. One is 
associated with the potential energy of raised liquid drops - 
this potential energy dissipates to friction as the drops fall. 
The second form of potential energy is associated with the 
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formation of a non-equilibrium pressure gradient, as the re- 
moval of vapor from the gas phase creates a pressure short- 
age of moist air aloft. This pressure gradient produces air 
movement. In the stationary case total frictional dissipation 
in the resulting circulation is balanced by the fraction of solar 
power spent on the work associated with evaporation. 

Evaporation is a surface-specific process. It is predomi- 
nantly anchored to the Earth's surface. In the stationary case, 
as long there is a supply of energy and the relative humidity 
is less than unity, evaporation is adding water vapor to the 
atmospheric column without changing its temperature. The 
rate of evaporation is affected by turbulent mixing and is usu- 
ally related to the horizontal wind speed at the surface. The 
global mean power of evaporation cannot exceed the power 
of solar radiation. 

In contrast, condensation is a volume-specific, rather than 
an area-specific, process that affects the entire atmospheric 
column. The primary cause of condensation is the cooling of 
air masses as the moist air ascends and its temperature drops. 
Provided there is enough water vapor in the ascending air, at 
a local and short-term scale condensation is not governed by 
solar power but by stored energy and can occur at an arbitrar- 
ily high rate dictated by the vertical velocity of the ascending 
flow, see p4] ). 

Any circulation pattern includes areas of lower pressure 
where air ascends, as well as higher pressure areas where it 
descends. Condensation rates are non-uniform across these 
areas - being greater in areas of ascent. Importantly, in such 
areas of ascent condensation involves water vapor that is lo- 
cally evaporated along with often substantial amounts of ad- 
ditional water vapor transported from elsewhere. Therefore, 
the mean rate of condensation in the ascending region of 
any circulation pattern is always higher than the local rate 
of evaporation. This inherent spatial non-uniformity of the 
condensation process determines horizontal pressure gradi- 
ents. 

Consider a large-scale stationary circulation where the re- 
gions of ascent and descent are of comparable size. A rel- 
evant example would be the annually averaged circulation 
between the Amazon river basin (the area of ascent) and 
the region of Atlantic ocean where the air returns from the 
Amazon to descend depleted of moisture. Assuming that 
the relative humidity at the surface, horizontal wind speed 
and solar power are approximately the same in the two re- 
gions, mean evaporation rates should be roughly similar as 
well (i.e., coincide at least in the order of magnitude). How- 
ever, the condensation (and precipitation) rates in the two re- 
gions will be consistently different. In accordance with the 
picture outlined above, the average precipitation rate P a in 
the area of ascent should be approximately double the av- 
erage value of regional evaporation rate E a . The pressure 
drop caused by condensation cannot be compensated by lo- 
cal evaporation to produce a net zero effect on air pressure. 
This is because in the region of ascent both the local wa- 
ter vapor evaporated from the forest canopy of the Amazon 



forest at a rate E a ~ Ed as well as imported water vapor 
evaporated from the ocean surface at a rate Ed precipitate, 
Pa = Ed + E a . This is confirmed by observations: precipi- 
tation in the Amazon river basin is approximately double the 
regional evaporation, P a ss 2E a (Marengo 2004 1. The dif- 



ference between regional rates of precipitation and evapora- 
tion on land, R = P a — E a ~ E a , is equal to regional runoff. 
Note that in the region of descent the runoff thus defined is 
negative and corresponds to the flux of water vapor that is ex- 
ported away from the region with the air flow. Where runoff 
is positive, it represent the flux of liquid water that leaves the 
region of ascent to the ocean. 

The fact that the climatological means of evaporation and 
precipitation are not commonly observed to be equal has 
been recognized in the literature (e.g., Wacker a nd Herbert] 
|2003j >, as has the fact that local mean precipitation values 
are consistently larger than those for evaporation (e.g., Tren 
|berthetal.|[200l] l. 

The inherent spatial non-uniformity of the condensation 
process explains why it is condensation that principally de- 
termines the pressure gradients associated with water vapor. 
So, while evaporation is adding vapor to the atmosphere 
and thus increasing local air pressure, while condensation 
in contrast decreases it, the evaporation process is signifi- 
cantly more even and uniform spatially than is condensation. 
Roughly speaking, in the considered example evaporation in- 
creases pressure near equally in the regions of ascent and 
descent, while condensation decreases pressure only in the 
region of ascent. Moreover, as discussed above, the rate at 
which the air pressure is decreased by condensation in the 
region of ascent is always higher than the rate at which lo- 
cal evaporation would increase air pressure. The difference 
between the two rates is particularly marked in heavily pre- 
cipitating systems like hurricanes, where precipitation rates 
associated with strong updrafts can exceed local evaporation 



rates by more than an order of magnitude (e.g., Trenberth and 
|Fasullo||2007] l. 

We have so far discussed the magnitude of pressure gradi- 
ents that are produced and maintained by condensation in the 
regions where the moist air ascends. This analysis is applica- 
ble to observed condensation processes that occur on differ- 
ent spatial scales, as we illustrated on the example of Hadley 
Cell. We emphasize that to determine where the ascending 
air flow and condensation will predominantly occur is a sep- 
arate physical problem. For example, why the updrafts are 
located over the Amazon and the downdrafts are located over 
the Atlantic ocean and not vice versa. Here regional evap- 
oration patterns play a crucial role. In Section 4.1 we have 
shown that constant relative humidity associated with surface 
evaporation, which ensures that dN v /dr = 0, is necessary 
for the condensation to take place. Using the definition of 
7 Q equation (37 » can be re-written as follows: 



d In 7 
Or 



w dj 
u dz 



(42) 
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This equation shows that the decrease of 7 with height and, 
hence, condensation is only possible when 7 grows in the 
horizontal direction, dlnj/dr > 0. Indeed, surface pressure 
is lower in the region of ascent. As the air moves towards the 
region of low pressure, it expands. In the absence of evapo- 
ration, this expansion would make the water vapor contained 
in the converging air unsaturated. Condensation at a given 
height would stop. 

Evaporation adds water vapor to the moving air to keep 
water vapor saturated and sustain condensation. The higher 
the rate of evaporation, the larger the ratio w/u at a given 



dj/dz and, hence, the larger the pressure gradient (37 1 that 



can be maintained between the regions of ascent and de- 
scent. A small, but persistent difference in mean evapora- 
tion AE < E between two adjacent regions, determines the 
predominant direction of the air flow. This explains the role 
of the high leaf area index of the natural forests in keeping 
evaporation higher than evaporation from the open water sur- 
face of the ocean, for the forests to become the regions of low 
pressure to draw moist air from the oceans and not vice versa 
(Mak arieva and Gors hkov, 2007). On the other hand, where 
the surface is relatively homogeneous with respect to evap- 
oration (e.g., the oceanic surface), the spatial and temporal 
localization of condensation events can be of random nature. 



5 Discussion: Condensation dynamics versus differen- 
tial heating in the generation of atmospheric circula- 
tion 

In Section 2 we argued that condensation cannot occur adi- 
abatically at constant volume but is always accompanied by 
a pressure drop in the local air volume where it occurs. We 
concluded that the statement that "the pressure drop by adi- 
abatic condensation is overcompensated by latent heat in- 
duced pressure rise of the air" (Poschl 2009 p. S 12437) was 
not correct. In Section 3 we quantified the pressure change 
produced by condensation as dependent on altitude in a col- 
umn in hydrostatic balance, to show that in such a column the 
pressure drops upon condensation everywhere in the lower 
atmosphere up to several kilometers altitude, Fig. lc. The 
estimated pressure drop at the surface increases exponen- 
tially with growing temperature and amounts to over 20 hPa 
at 300 K, Fig. lb. 

In Section 4 we discussed the implications of the 
condensation-induced pressure drop for atmospheric dynam- 
ics. We calculated the horizontal pressure gradients produced 
by condensation and the efficiency of the atmosphere as a dy- 
namic machine driven by condensation. Our aim throughout 
has been to persuade the reader that these implications are 
significant in numerical terms and deserve a serious discus- 
sion and further analysis. We will now conclude our consid- 
eration by discussing the condensation-induced dynamics at 
the background of differential heating, a physical mechanism 



that, in contrast to condensation, has received much attention 
as an air circulation driver. 

Atmospheric circulation is only maintained if, in agree- 
ment with the energy conservation law, there is a pressure 
gradient to accelerate the air masses and sustain the existing 
kinetic energy of air motion against dissipative losses. For 
centuries, starting from the works of Hadley and his prede- 
cessors, the air pressure gradient has been qualitatively asso- 
ciated with the differential heating of the Earth's surface and 
the Archimedes force (buoyancy) which makes the warm and 
light air rise, and the cold and heavy air sink. This idea can 
be illustrated by Fig. lc, where the warmer atmospheric col- 
umn appears to have higher air pressure at some heights than 
the colder column. In the conventional paradigm, this is ex- 
pected to cause air divergence aloft away from the warmer 
column, which, in its turn, will cause a drop of air pressure at 
the surface and the resulting surface flow from the cold to the 
warm areas. Despite the physics of this differential heating 
effect being straightforward in qualitative terms, the quanti- 
tative problem of predicting observed wind velocities from 
the fundamental physical parameters has posed enduring dif- 
ficulties. Slightly more than a decade before the first signif- 



icant efforts in computer climate modelling, Brunt ( 19441 as 



cited by Lewis ( 1998] > wrote: 

"It has been pointed out by many writers that it is impossible 
to derive a theory of the general circulation based on the known 
value of the solar constant, the constitution of the atmosphere, and 
the distribution of land and sea . . . It is only possible to begin 
by assuming the known temperature distribution, then deriving the 
corresponding pressure distribution, and finally the corresponding 
wind circulation". 

Brunt's difficulty relates to the realization that pressure 
differences associated with atmospheric temperature gradi- 
ents cannot be fully transformed into kinetic energy. Some 
energy is lost to thermal conductivity without generating me- 
chanical work. This fraction could not be easily estimated by 
theory in his era - and thus it has remained to the present. 
The development of computers and appearance of rich satel- 
lite observations have facilitated empirical parameterizations 
to replicate circulation in numerical models. However, while 
these models provide reasonable replication of the quantita- 
tive features of the general circulation they do not constitute a 
quantitative physical proof that the the observed circulation 
is driven by pressure gradients associated with differential 



heating. As Lorenz (1967 p. 48) emphasized, although "it 
is sometimes possible to evaluate the long-term influence of 
each process affecting some feature of the circulation by re- 
course to the observational data", such knowledge "will not 
by itself constitute an explanation of the circulation, since it 
will not reveal why each process assumes the value which it 
does". 

In comparison to temperature-associated pressure differ- 
ence, the pressure difference associated with water vapor 
removal from the gas phase can develop over a surface of 
uniform temperature. In addition, this pressure difference 
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is physically anchored to the lower atmosphere. Unlike the 
temperature-related pressure difference, it does not demand 
the existence of some downward transport of the pressure 
gradient from the upper to the lower atmosphere (i.e., the 
divergence aloft from the warmer to the colder column as 
discussed above) to explain the appearance of low altitude 
pressure gradients and the generation of surface winds. 

Furthermore, as the condensation-related pressure differ- 
ence Sp s is not associated with a temperature difference, 
the potential energy stored in the pressure gradient can be 
nearly fully converted to the kinetic energy of air masses 
in the lower atmosphere without losses to heat conductiv- 
ity. This fundamental difference between the two mecha- 
nisms of pressure difference generation can be traced in hur- 
ricanes. Within the hurricane there is a marked pressure gra- 
dient at the surface. This difference is quantitatively account- 



able by the condensation process (Makarieva and Gorshkov 
2009b R In the meantime, the possible temperature differ- 



ence in the upper atmosphere that might have been caused 
by the difference in moist versus dry lapse rates between the 
regions of ascent and descent is cancelled by the strong hor- 
izontal mixing (Montgomery et al. 2006| l. Above approx- 
imately 1.5 km the atmosphere within and outside the hur- 
ricane is approximately isothermal in the horizontal direc- 
tion ( [Montgomery et aL| |2006| Fig. 4). Therefore, while the 
temperature-associated pressure difference above height z c , 
Fig. lc, is not realized in the atmosphere, the condensation- 
associated pressure difference below height z c apparently is. 

Some hints on the relative strengths of the circulation 
driven by differential heating compared to condensation- 
induced circulation can be gained from evaluating wind ve- 
locities in those real processes that develop in the lower at- 
mosphere without condensation. These are represented by 
dry (precipitation-free) breezes (such as diurnal wind pat- 
terns driven by the differential heating of land versus sea sur- 
faces) and dust devils. While both demand very large temper- 
ature gradients (vertical or horizontal) to arise as compared 
to the global mean values, both circulation types are of com- 
paratively low intensity and have negligible significance to 
the global circulation. For example, dust devils do not in- 
volve precipitation and are typically characterized by wind 
velocities of several meters per second ( |Sinclair||1973| l. The 
other type of similarly compact rotating vortexes - tornadoes 

- that are always accompanied by phase transitions of water 

- develop wind velocities that are at least an order of magni- 
tude higher ( |Wurman et al. 1996k More refined analyses of 
Hadley circulation (Hel d and Hou|[T980| l point towards the 
same conclusion: theoretically described Hadley cells driven 
by differential heating appear to be one order of magnitude 
weaker than the observed circulation ( Held and Hou| |1980[ 
|Schneider| [20061 ), see also ( |Caballero efaTl|2008) . While the 
theoretical description of the general atmospheric circulation 
remains unresolved, condensation-induced dynamics offers a 
possible solution (as shown in Section 4.1). 

Our approach and theory have other significant impli- 



cations. Some have been documented in previous papers, 
for example with regard to the development of hurricanes 
( Makarieva and Gorshkov 2009a|b i and the significance of 
vegetation and terrestrial evaporation fluxes in determining 
large scale continental weather patterns (Makarieva et al. 



2006 ; Makarieva and Gorshkov] |2007[[Sheirand Murdiyarso 



2009| |Makarieva et al. 2009). Other implications are likely 



to be important in predicting the global and local nature of 
climate change - a subject of considerable concern and de- 
bate at the present time ( |Pielke et al.| |2009| |Schiermeier| 
|2OT0l 

In summary, although the formation of air pressure gradi- 
ents via condensation has not received detailed fundamental 
consideration in climatological and meteorological sciences, 
here we have argued that this lack of attention has been un- 
deserved. Condensation-induced dynamics emerges as a new 
field of investigations that can significantly enrich our under- 
standing of atmospheric processes and climate change. We 
very much hope that our present account will provide a spur 
for further investigations both theoretical and empirical into 
these important, but as yet imperfectly characterized, phe- 
nomena. 
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